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1. Introduction 

The LHC searches for jets +$t events, as well 
as searches for heavy Higgs bosons H/A and the 
rare decay B s — > , are placing ever-stronger 

constraints on supersymmetric models, in partic- 
ular those in which R-parity is conserved and the 
lightest neutralino Xi is the stable lightest su- 
persymmetric particle and provides the cosmo- 
logical dark matter. There have been a number 
of analyses of the constraints imposed by AT- 
LAS and CMS using ~ 1/fb at 7 TeV [Iffi and 
other data on the parameter spaces of simpli- 
fied versions of the minimal supersymmetric ex- 
tension of the Standard Model (MSSM) [|, par- 
ticularly the constrained MSSM (CMSSM) [t] in 
which the soft supersymmetry-breaking parame- 
ters tuq, toi/2 and Aq are universal at the GUT 
scale, and some also of the LHC constraints on 
models in which the soft supersymmetry-breaking 
contributions to the masses of the Higgs multi- 
plets are non-universal but equal to each other, 
the NUHM1 models [fjj. There have also been 
analyses of the CMSSM and NUHM1 of the 
prospective implications of a Higgs boson weigh- 
ing ~ 125 GeV [2|[7}fT0], and analyses within 



the CMSSM and NUHM1 [8f|ll| that incorporate 
LHC searches for jets + 1/Jt events with 
of LHC data at 7 TeV fl2l[l3 



5/fb 

Here we extend 
our previous frequentist analyses of the CMSSM 
and NUHM1 [TJ[2j[l4}|20] to include the ATLAS 
jets +lpT results with ~ 5/fb [l2j of data as well 
as Mh ~ 125 GeV and recent data on BR(B S — ► 
M+aO from the ATLAS (2l], CDF [22], CMS [23 
and LHCb 24 Collaborations [25], a nd the recent 



fi+fi 



2G 



We also im- 



CMS exclusion of H/A 
plement the constraint on the spin-independent 
dark matter scattering cross section imposed by 
225 live days of XENON100 data [27]. We give 



a critical discussion of the extent to which the 
CMSSM and NUHM1 can be said to be dis- 
favoured by the available data from the LHC and 
elsewhere. 

The results of experimental searches for jets 
+ J/jt events are typically presented within the 
framework of the CMSSM for some fixed A and 
tan j3. The applicability of these analyses to other 
Aq and tan/3 values, as well as to constraining 



the NUHM1, requires some study and justifica- 
tion. Within the CMSSM, one must verify that 
the constraints used are indeed independent of Aq 
and tan/3, as often stated in the experimental pa- 
pers on searches for jets + T/t events. Similarly, 
within the NUHM1 one must verify whether, 
for any specific set of values of too, TOi/2j Aq 
and tan j3, the sensitivities of ATLAS and CMS 
searches for jets +1/t events might depend on the 
degree of non- universality in the NUHM1. A sec- 
ond issue arising in the NUHM1 is that the range 
of mo that is consistent with the Xi LSP require- 
ment depends on the degree of non-universality. 
In the CMSSM, for any given value of m\i<i there 
is a lower limit on mo that is, in general, vio- 
lated in the NUHM1. Since the LHC experiments 
generally quote results only for the region of the 
(too,TOi/2) plane allowed within the CMSSM, a 
dedicated study is needed to estimate the correct 
extrapolation of LHC results outside this region. 

In order to explore the sensitivity to more gen- 
eral CMSSM and NUHM parameters of the pub- 
lished LHC searches for supersymmetry within 
specific CMSSM models, we use a version of the 
Delphes [28] generic simulation package with a 
'card' that emulates the performance of the AT- 
LAS detector. As described below, in the specific 
case of the 7 TeV ATLAS 5/fb jets + $ T anal- 
ysis, Delphes reproduces reasonably accurately 
the quoted 95% CL limits on CMSSM models 
with Aq = and tan/3 = 10, as well as other con- 
fidence levels in the CMSSM (too,TOi/2) plane. 
These are well described by scaling of event num- 
bers oc 1/M 4 , where M. = ^jm\ + to^ 2 , as as- 
sumed in [l]. We also use the Delphes simulation 
to confirm that the ATLAS 5/fb jets + 1/ T con- 
straint in the {mo, mi/2) plane of the CMSSM is 
relatively insensitive to tan/3 and Aq. 

Having validated the application of the 
Delphes simulation to the CMSSM, we then 
study whether the results of 12 1 are sensitive 
to the degree of Higgs non-universality in the 
NUHM1. Although the supersymmetric decay 
cascades to which this analysis is most sensitive 
are independent of the details of the supersym- 
metric Higgs sector, they do in principle have 
some sensitivity to the Higgs mixing parameter /z, 
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which is in general different in the NUHMf from 
its value in the CMSSM. We find that in prac- 
tice the jets + ]£t constraint is essentially inde- 
pendent of the degree of non-universality in the 
NUHM1. We then use Delphes to estimate the 
sensitivity of the ATLAS analysis to the region 
with mo smaller than is allowed in the CMSSM, 
finding that the ATLAS jets + I^t constraint is 
essentially independent of too in this range. 

Following these validations of our implemen- 
tation of Delphes in the NUHM1, we use the 
Markov Chain Monte Carlo (MCMC) technique 
to sample the CMSSM and NUHM1 parameter 
spaces with ~ 95 x 10 6 and ~ 221 x 10 6 points, 
respectively. Given the constraints from Mh and 
the LHCi/fb data, we find that the ATLAS 5/fb 
jets + IfjT constraint by itself has a relatively lim- 
ited impact on the 68 and 95% CL regions in the 
(TO ,m 1/2 ) planes of the CMSSM and NUHM1, 
whereas the new BR(£> S — > constraint has 

larger effects in both these and the (Ma, tan (3) 
planes, and the new XENON100 constraint also 
impacts significantly both the (toq,TOx/2) and 
(M A ,tan/3) planes of the CMSSM and NUHM1. 
After the inclusion of these constraints, we find 
that in the favoured regions in the CMSSM and 
NUHM1 in mo, miw and tan/3 the global x 2 
function varies little, giving less relevance to the 
best-fit values of these parameters. Nevertheless, 
we note that the best-fit values of to and TOi/ 2 
are decreased in the CMSSM, and the best-fit val- 
ues of tan /3 are decreased in both the CMSSM 
and the NUHM1. Within the CMSSM, the global 
X increases from 31.0 to 32.8 and the p-value 
is correspondingly reduced from 12% to 8.5%. 
Within the NUHM1, the value of tan/3 rises to a 
value similar to that in the CMSSM, while the to- 
tal x 2 rises from 28.9 to 31.3, reducing the p- value 
from 15% to 9.1%. The perilously low p-values of 
the CMSSM and NUHM1 reflect the growth in 
tension between LHC constraints and low-energy 
measurements such as (g 
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2. Analysis Procedure 



We follow closely the procedure described in [T] . 
As in [l] , the sampling is based on the Metropolis- 
Hastings algorithm, with a multi-dimensional 



Gaussian distribution as proposal density. The 
width of this distribution is adjusted during the 
sampling, so as to keep the MCMC acceptance 
rate between 20% and 40% in order to ensure ef- 
ficient sampling. We emphasize that we do not 
make use of the sampling density to infer the un- 
derlying probability distribution. Rather, we use 
the MCMC method to construct a global like- 
lihood function that receives contributions from 
the usual full set of electroweak precision observ- 
ables including (g — 2) M [29], measurements of 
BR(& -> s 7 ), BR(B U -> rz/ T ), BR(B S -> m + aO, 
B-meson mixing and other flavour-physics observ- 
ables [25pi] , the cosmological Xi density and the 
XENON100 direct search for dark matter scatter- 
as well the LHC searches for the Higgs 
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boson and supersymmetric signals. 

Our treatment of the non-LHC constraints is 
similar to that in 1|, with the exceptions that 
we use updated values of Mw — 80.385 ± 
0.015 GeV 
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and m t = 173.2 ± 0.9 GeV 
Our treatment of BR(_B S — > is based on a 

compilation of the recent measurements by AT- 
CDF [22" 
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CMS [23] and LHCb [24][25]. 
As already mentioned, we incorporate the public 
results of the search for jets + ]£t events with- 
out leptons using ~ 5/fb of LHC data at 7 TeV 
analyzed by ATLAS [12], which has greater sen- 
sitivities to the models discussed here than do 
searches with leptons and 6-quarks in the final 
state as well as searches with less integrated lu- 
minosity. We also include the CMS constraint 
on the heavier MSSM Higgs bosons, H/A [26] , 
which has greater sensitivity than the correspond- 
ing ATLAS search [35]. 

Concerning the Higgs boson, the situation has 
evolved rapidly since our previous analysis J2] 



and the data reported earlier this year 39 40 



with the recently-announced discovery by the 
CMS 
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con- 



and ATLAS Collaborations 
sistent with evidence from the CDF and DO Col- 
laborations 
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of a new particle that is a very 
strong candidate to be a Higgs boson resembling 
that in the Standard Model. In the absence of 
an official combination of these results, we model 
them by assuming, as in (2), that 



M h = 125 ± 1 (exp.) ± 1.5 (theo.) GeV , 



(1) 
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and we incorporate this new constraint using the 
same 'afterburner' approach as in [2]. 

The contributions of these observables to the 
global likelihood function are calculated within 
the MasterCode framework 



passed through the Delphes 28 simulation code 
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This incor- 
porates a code for the electroweak observables 
based on [42] as well as the SoftSUSY [43 



FeynHiggs [44], SuFla [45], and MicrOMEGAs [46 

codes, interfaced via the SUSY Les Houches 
Accord 



47 



48 and 



We use the Super I so 
SSARD |49]_c odes as comparisons for B-physics ob- 
servables]^] and for the cold dark matter density, 
using also DarkSUSY [50] as a comparison for the 
latter. We use SSARD for the spin-independent 
dark matter scattering cross section, comparing 
with MicrOMEGAs and DarkSUSY. 

3. Validation and Extension of Jets + $t 
Constraints 

ATLAS and CMS results on SUSY searches for 
jets +E T ar e typically shown as 95% CL exclusion 
bounds in the CMSSM {m ,m 1 / 2 ) plane for fixed 
tan/3 and A , e.g., tan/3 = 10 and A = 12 
53 . We first validate for these values of tan (3 and 
Aq our implementation of the ATLAS constraint 



2* 



using PYTHIA ;54j and the generic Delphes 
simulation code with an ATLAS detector 'card' 
Next we extend it to other values of the CMSSM 
parameters, and then to the NUHM1 

Validation in the (mo, mi/2) plane 

As a first step in the validation procedure, 
we used SoftSUSY [43] to generate the spectra 
for various points in the (mojm^) plane of the 
CMSSM for tan/3 = 10 and A = 0. We then 



used PYTHIA 54 to generate 50,000 events for 
each of these points, which were subsequently 



We do not inclu de the isospin asymmetry in B — > 
K^*) u + u— d ecays |5l| or the measurement of BR(B — > 
D^tu) 52] in our analysis, in view of the present exper- 
imental and theoretical (long-distance) uncertainties. 
2 We have made a similar validation analysis for the CMS 
olt search for jets + $t events using ~ 1/fb of data at 
7 TeV [53| . We do not discuss this validation in detail, as 
it does not contribute to the the likelihood function ana- 
lyzed here, but it does validate a posteriori our previous 
treatments [l][2] of the CMS analysis. It also indicates 
that the CMS sensitivity with 5/fb of data at 7 TeV [lj] is 
similar to the ATLAS 5/fb jets + $t data discussed here. 



using the ATLAS detector 'card'. We then used 
information from 112] to calculate the efficien- 
cies for observing events that would survive the 
cuts used in the ATLAS jets + J$t analysis us- 
ing 5/fb of data at 7 TeV. Taking into account 
the observed number of events, together with the 
information in Table 3 of [12] , we then calcu- 
lated the confidence levels of these points, us- 
ing the same prescription as ATLAS. We recall 
that the ATLAS analysis combined searches us- 
ing several different channels A, A',B,C, D and 
E, which, depending on the channel, may have 
'tight', 'medium' and 'loose' event selections. Fol- 
lowing the procedure used in the ATLAS analysis, 
for each set of CMSSM parameters we used the 
limit set by the most sensitive of these various 
channel selections. 

The left panel of Fig. [T] displays how closely 
this work flow reproduced the ATLAS CL for 
34 test points along its 95% CL exclusion line, 
sampling the range m G (200, 3500) GeV in 
100 GeV steps. Our results indicate a mean dif- 
ference CL Delphes - CL ATL as = (-2.8 ± 3.6)%, 
so we conclude that our simulation using PYTHIA 
and Delphes reproduces quite accurately the 95% 
CL exclusion contour in the (mo,mi/2) plane for 
tan /3 = 10 and A = reported by ATLAS. 

Since it is not computationally feasible to use 
PYTHIA and Delphes in this way for every point 
generated in our MCMC analysis, we make use of 
a simple analytic approximation to their results 
in order to extrapolate their results to other re- 
gions of the (wo, 7714/2) plane. The right panel 
of Fig. [T] compares the CL values estimated us- 
ing this Delphes implementation, for a number 
of CMSSM points with tan ,3 = 10, A Q = and 
varying mo and mi/2, with the CL calculated as- 
suming, as in [l], that the numbers of events for 
different points in the CMSSM {m a ,m 1 / 2 ) plane 
scale as 1/A4 4 (solid line). We see that the agree- 
ment is quite good along rays in the (mo 1 rn 1 / 2 ) 
plane with mi/2/mo = 1/3 (blue circles), 1 (red 
squares) and 3 (green triangles), indicating that 
this scaling law is an adequate approximation to 
the true sensitivity of the ATLAS 5/fb jets +E T 
analysis. 
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Figure 1. Left panel: A histogram of the differences between the confidence levels (CLs) we calculate 
using PYTHIA and Delphes from the values quoted by ATLAS for 34 points along the 95% exclusion 
contour in the (7710,^1/2) plane of the CMSSM with tan/3 = 10 and Ao = 0. The mean difference is 
C^Deiphes - CL atlas = (—2.8 ±3.6)%. Right panel: The points are the confidence levels for a selection of 
CMSSM points with tan /3 = 10 and Aq = and various values of R = A4/A4q, where A4 = ^Jm^ + m^ 2 

and Mq is the value of M. at the point on the LHC jets + 1£t 95% exclusion line with the same ratio 
mo/m-L/2- The red squares (blue circles, green triangles) were chosen along rays in the (mo, mm) plane 
with mo/mi/2 = 1/3,1,3, respectively. These points compare well with the solid line, which is the CL 
calculated assuming that the number of signal events surviving the LHC cuts scales as 1/Ai 4 . 



Extension to other tan fj and A values 

The next step is to use our Delphes imple- 
mentation to evaluate the sensitivity of the AT- 
LAS analysis to other values of tan (3 and Ao . 
To this end, for each of the 34 points along 
the 95% CL contour in the (mo,mi/ 2 ) plane for 
tan/3 = 10 and Aq = mentioned previously, we 
used PYTHIA and Delphes to estimate the confi- 
dence levels of these points as the other parame- 
ters are varied, as shown in Fig. [2] The left panel 
exhibits the variation of the confidence level with 
tan /3 for fixed Aq = 0. We see that the CL is quite 
stable, with its central value ranging between 95 
and 90% and its 68% CL range lying between 
99% and 84% and always including the 95% CL 



found by ATLAS for tan/3 = 10 and A = Q 
The right panel of Fig. [2] displays a similar anal- 
ysis for varying Ao/rug € [—3, 3] with tan/3 fixed 
at 10. We see again that the CL is quite stable, 
with central values between 95 and 92% and the 
68% CL ranges lying between 100% and 88% and 
again always including the 95% CL found by AT- 
LAS for tan /3 = 10 and Ao = 0. These results are 
consistent with the statement by ATLAS [12] that 
their (mo, miw) exclusion contours are largely in- 
dependent of tan /3 and A , as we have assumed 
previously [1 j , and also provide consistency checks 
on our Delphes implementation. 

3 The different sizes of the error bars, here and in subse- 
quent validation plots, are due to PYTHIA failing for varying 
numbers of points. 
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Figure 2. Calculations using Delphes of the average and spread in exclusion confidence levels of 34 
CMS SM points selected to have the same (too,toi/2) values as points on the 95% ATLAS 5/fb jets +Ipr 
exclusion line for tan /3 = 10 and Aq = 0, but with different values o/tan/3 (left panel) and Aq (right 
panel). The CL is almost independent of these parameters, with the error bar always embracing 95%. 
The band represents the Ha range of the fit shown in the left panel of Fig. [7} 



Extensions to Non-Universal Higgs Masses 

Following these validations, as a third step 
we now use our Delphes implementation to ex- 
tend the interpretation of the ATLAS 5/fb jets + 
Ijjf analysis to models with non-universal Higgs 
masses. We recall that in the CMSSM, the elec- 
troweak vacuum conditions can be used to fix the 
Higgs mixing parameter /i and Ma for any given 
choice of the parameters mo,mi/2, A and tan/3. 
Introducing the single extra degree of freedom 
m H / m o = m # / m o m the NUHM1, one can 
treat either \x or Ma as a free parameter, whereas 
these quantities may be treated as two free pa- 
rameters in the NUHM2, in which m 2 H /toq and 
m 2 Hd /rn^ are treated as two independent param- 
eters [55) . The effects of varying fi and Ma on 
the gluino mass and the spectra of squarks (the 
main sparticles expected to be produced at the 



LHC) are of secondary importance [^] However, 
they may also affect the branching fractions of the 
cascade decays to which the LHC experiments are 
sensitive, and so could in principle have some ef- 
fect on the sensitivity of the ATLAS 5/fb jets + 
$t search in the (m ,mi/ 2 ) plane. 

For this reason, we have taken the same 
34 points along the 95% CL contour in the 
{m ,m 1 / 2 ) for tan/3 = 10 and A = studied 
above, and used PYTHIA and Delphes to estimate 
how the confidence levels of these points vary as 
the Higgs non-universality parameters m 2 H /m§ 
and w^/toq are varied. The red squares in 
the left panel of Fig. [3] show how the CL de- 
pends on the common non-universality param- 
eter m| o /mQ = rrL 2 Hd /mQ in the NUHM1. As 
before, we see that the CL does not depart sig- 



4 Except that varying affects the masses and mixings of 
third-generation squarks, but the dominant LHC jets + 
IpT searches are less sensitive to these. 
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nificantly from the CMSSM value of 95%. Also 
shown in the left panel of Fig. [3] (as blue circles) 
are PYTHIA and Delphes estimates of the con- 
fidence levels of points with the same values of 
Wo, mi/2, tan/3 = 10 and Aq = but with Higgs 
non-universality parameters of opposite signs for 
the two Higgs doublets: m% /m 2 , = —rn 2 Hd /rn^. 
The right panel of Fig. [3] shows how the CL in 



the NUHM2 varies as a function of ; 



H, 



for fixed 



m 2 H = (red squares) and as a function of m 2 H 
for fixed fn\ d = (blue circles). We again see 
that the CL is similar to the 95% found in the 
CMSSM for the same values of m , roirt, tan j3 
Hi and A Q = 0. 

Extension to m 2 , < 

To complete the validation for the NUHM1, 
and in preparation for future analyses in the 
NUHM2, we have also investigated the applicabil- 
ity of the ATLAS 5/fb jets +f T results to models 
with smaller values of m\ than in the CMSSM, 
including negative values. We have evaluated the 
CL along a line of points with tan/3 = 10, Aq = 
0,m 1/2 = 630 GeV and -(150 GeV) 2 < m 2 < 
+(150 GeV) 2 , i.e., M = sign(m^) x y/\m%\ G 
(-150, +150) GeV, as shown in Fig. |IJ We see 
that the CL remains constant, within errors, over 
this range of Mo, and is quite consistent with the 
95% CL reported for the ATLAS 5/fb jets +$ T 
analysis at m = 100 GeV. 

4. Combination of BK(B S /i + /i~) Con- 
straints 

Since our previous analysis [2] , new constraints 
on BR(B S — > have been published by AT- 

LAS " 



121], 



24 



CDF 122], CMS [23] and LHCb 
Collaborations [25] , Their measurements are to 
be compared with the theoretical (TH) value 
in the SM, which we take to be BR(_B S — > 
M + /Osm,th = (3.2 ± 0.2) x 10~ 9 [57]. This 
comparison is not direct, because experiments 
measure the time-averaged (TA) branching ra- 
tio, which differs from the TH value because of 
the difference between the lifetimes of the heav- 
ier and lighter B s mesons [58] [59] . In general, 
BR(B S -> m + aOth = [(1 -H#7(l +A A ry s )] x 
BR(_B S — > /i + /i~)ta, where ^Ar = +1 in the 



SM and y s = 0.088 ± 0.014. Thus we com- 
pare the measurements to the SM prediction for 
the time-averaged branching ratio: BR(_B S — > 
/i + /OsM.TA = BR(_B S ^>- /x+^") S m.th/(1 
y s ) = (3.5 + 0.2) x 10- 
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We use for this comparison an unofficial com- 
bination of the ATLAS, CDF, CMS and LHCb 
data 21 ■ 25] , from which we compute a \ 2 penalty 
for each possible (calculated) value of BR(i3;5 — > 
/i + /i _ ). In Fig. [5] we show the probability 
distribution (left) and the log-likelihood func- 
tion for the ratio BR(£? S — > ^ + /i")/BR(B s — > 
/x + (U _ )sm,ta, folding together the theoretical and 
experimental errors. The log-likelihood function 
is used as the A% 2 penalty for the values of this 
observable as they are calculated in the CMSSM 
and NUHM1 

5. The Constraint on Spin-Independent 
Dark Matter Scattering from 225 Live 
Days of XENON100 Data 

The XENON100 Collaboration has recently 
published the result of a search for spin- 
independent dark matter scattering based on 225 
live days of data (27 . They report two events 
in a fiducial region where 1.0 ± 0.2 background 
events were expected, i.e., no evidence for any 



si 



< 



excess, quote a 90% CL upper limit of 
2.0 x 10 -45 cm 2 for a dark matter particle mass 
of 55 GeV, and provide a 90% CL exclusion line 
for other masses. 



Following the procedure used in 20 for the pre- 
vious XENON100 constraint, we use a CL S cal- 
culator, now implemented in our Delphes frame- 
work, to model the contribution to the global x 2 
function of this new XENON100 constraint from 
the number of signal events as a Gaussian with 
mean /i = 1.0 and standard deviation a — 2.7 
events. This corresponds to the mean (2 - 1 = 
1 event) given in 



27 and the upper limit of 5.1 



signal events at the 90% CL given by the CL S 
calculator, which we identify with the published 
XENON100 90% exclusion line. For models pre- 
dicting values of a^ 1 above or below the 90% 
CL line in [27], we rescale linearly the number 



5 We note in passing that one expects _4at = +1 also in 
these models. 
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Figure 3. Calculations using Delphes of the average and spread in exclusion confidence levels of 34 
CMS SM points selected to have the same (too, m i/2) values as points on the 95% LHC jets +$t exclusion 
line for tan/3 = 10 and Aq = 0, as in Fig. Q but (left panel) varying rn 2 Hd /m'( i = np? H ^lrn\ in the NUHM1 
(red squares) and vtv 2 Hd jm^ = —m 2 H /toq in the NUHM2 (blue circles), and (right panel) varying m 2 Hd 
for fixed rn 2 H ^ = (red squares) and varying ra 2 H ^ for fixed rn 2 Hd — (blue circles) in the NUHM2. In 
both panels, the bands represent the Ha range of the fit shown in the left panel of Fig. [7} 



of events and estimate the A% 2 contribution of 
XENON100 using the Gaussian model described 
above. 

We use SSARD 491 to calculate ojj 1 , assuming a 
nominal value of the 7r-N a term E^tv = 50 MeV, 
and assigning an asymmetric error MeV. This 
accounts for the fact that E^at = ctq ~ 36 MeV 
if (iV|ss|iV) = 0, while allowing for significantly 
larger values of E^at as suggested by some anal- 
yses [60] . 

The XENON100 analysis assumes that the lo- 
cal dark matter density is 0.3 GeV/cm 3 , and 
makes supplementary assumptions on the dark 
matter velocity distribution. We have explored 
the potential implications of the uncertainties in 
these quantities by convoluting with our above- 
mentioned model of the XENON100 constraint 
on CTp 1 a fractional uncertainty of a factor \/2 
at the 68% CL. We find that it has negligible im- 



pact on the results shown below, which are for 
the nominal dark matter distribution parameters 
assumed by XENON100. 



6. The CMSSM and 



ATLASs/fb, 
XENON100 



BR(B S 



NUHM1 after 
') and 



As already mentioned, in addition to the 
Delphes implementation of the ATLAS 5/fb jets 
+ IpT constraint, the implementation of the com- 
bined BK(B S —> constraints described 
above, and the new CMS constraint on H/A — > 
t + t~ 26 , we include in our global fit the con- 
straint M h = 125 ± 1.0 ± 1.5 GeV and the other 
constraints included in [2|, with updated values 
of the W boson mass [33] and the top quark 
mass |34l, see above. 
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Figure 4. Calculations using Delphes of the exclusion confidence levels of NUHM1 points selected to have 
tan f3 = 10, Aq = 0,m 1/2 = (630 ± 25) GeV and varying M = sign{ml) x y/\m$[ £ (-150, +150) GeV. 
The band represents the ±lcr range of the fit shown in the left panel of Fig. [7J 
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Figure 5. The probability distribution (left panel) and the log-likelihood function (right panel) found 
in an unofficial combination of the constraints on HH(B S — > fi + fi~) from ATLAS, CDF, CMS and 
LHCb flip^ normalized to BR(B S -> m + /-Osm,ta = BR(B S -> M + M )sm,th/(1 - y s ) = (3.5 ± 0.2) x 
10~ 9 , as described in the text. 
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The (mo,mi/ 2 ) planes 

Fig. ^displays the (mo, 7774/2) planes for the 
CMSSM (left panel) and the NUHM1 (right 
panel). The results of the new global fits are 
shown as solid lines and filled stars, and previous 
fits based on ~ 1/fb of LHC data are indicated 
by dashed lines and open stars. The blue and red 
lines denote levels A\ 2 = 2.30 and 5.99, which 
we interpret as the 68% and 95% CLs. 

We note the existence in the CMSSM of a 
grouping of points in the focus-point region at 
large tuq and small 777,1/2 that were previously al- 
lowed at the 95% CL but are now disallowed [f] 
This exclusion is due to the new XENON100 
constraint [27]. The main effect of the ATLAS 
5/fb jets + Ipr constraint is to shave off a slice 
of points at low 777,1/2 and moderate mo. We 
also note the appearance in the NUHM1 sample 
of a grouping of points at large 777,1/2 that have 
smaller mo than is possible in the CMSSM, and 
are allowed at the 95% CL. The relic density in 
this region is brought into the WMAP range by 
multiple neutralino and chargino coannihilation 
processes. Points are found in this region using 
each of the codes MicrOMEGAs [46], DarkSUSY [50], 
Superlso [48] and SSARD [49], though they yield 
somewhat different values of the relic density. As 
mentioned previously, here we use MicrOMEGAs as 
our default. 

The values of x 2 in the most-favoured regions 
of the CMSSM and NUHM1 are relatively shal- 
low functions of mo, 7774/2 and tan/3, so that lit- 
tle importance can be attributed to the positions 
of the best-fit points or to their movements rela- 
tive to those in [2]. Nevertheless, we note that, 
as summarized in Table [T] the nominal best-fit 
point in the CMSSM moves to smaller mo and 
mi/2, whereas that in the NUHM1 moves little 
in this plane. The best-fit points have % 2 = 32.8 
in the CMSSM and 31.3 in the NUHM1, and the 
p-values of these models are reduced to 8.5 and 
9.1%, respectively. 

We see in Fig. [6] that the 95% CL region in 
the CMSSM extends to m ~ 4000 GeV and 
mi/2 ~ 2400 GeV. The corresponding region in 



6 Tliis could also be inferred from Figs. 7 and 9 of the paper 
by JE and KAO in [7|. 



the NUHM1 extends to smaller mo but larger 
TO1/2. We also note that the CMSSM fit fea- 
tures two disconnected 68% CL 'islands', the one 
at lower mo and 777,1/2 corresponding to the stau 
coannihilation region, and that at larger mo and 
7714/2 corresponding to the rapid-annihilation fun- 
nel region. We note, however, that our sam- 
pling also includes a few intermediate points 
that also have A% 2 < 2.30, reflecting the rela- 
tive flatness of the global x 2 function along the 
band between (mo, 7711/2) ~ (300,900) GeV and 
- (1100, 1900) GeV. Typically, these points sit in 
rapid-annihilation funnels, and it is possible that 
a more complete sampling might reveal a lower 
'isthmus' connecting these islands, as happens in 
the NUHM1. 

We report in Tablefllthe values of mo, mi/2, Aq 
and tan (3 for the best fits in both the CMSSM 'is- 
lands', as well as for the best fit in the NUHM1, 
which is at relatively low masses, and a local min- 
imum of the NUHM1 x 2 function at high masses. 
We note that the best fit in the low-mass CMSSM 
'island' is similar to the best fit in the NUHM1, 
and that both have smaller tan/3 than in the pre- 
vious best fits to the LHCi/fb data set in these 
models. As we discuss later, this effect is due 
to the new BK(B S —> /i + /j~) constraint. On the 
other hand, the best fit in the high-mass CMSSM 
'island' is more similar to the previous best fit to 
the LHCi/fb data set. In general terms, we see 
that the area of the NUHM1 (mo,m 1 / 2 ) plane 
that is allowed at the 68 and 95% CL is larger 
than in the (mo, 777,1/2) plane of the CMSSM. This 
reflects the fact that the dependence of the global 
X 2 function on these variables is weaker in the 
NUHM1 than in the CMSSM, reflecting the ad- 
ditional freedom offered by the degree of Higgs- 
mass non-universality. 

In particular, the 68% CL region in the 
NUHM1 is connected, and extends to larger m 
than in the CMSSM. This extension reflects the 
fact that in the NUHM1 a rapid-annihilation fun- 
nel region may appear at smaller 7774/2 and tan/3 
than in the CMSSM for the same value of mo, 
thanks to the variation in M4 that is possible in 
the NUHM1, allowing Ma ~ 2m^o in different re- 
gions of the (777,0,7774/2) plane. It is this freedom 
in the location of the rapid-annihilation funnel, a 
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Figure 6. The (m-Oi m i/2) planes in the CMSSM (left panel) and the NUHM1 (right panel) including 
the ATLAS 5/fb jets + $ T constraint a combination of the ATLAS J2^, CDF (M CMS {If 

and LHCb constraints on BR(_B S — > /i + /i~) '251 and the recent XENON100 result \27lj, assuming 
Mh = 125 ± 1 (exp.) ± 1.5 (theo.) GeV. The results of the current fits are indicated by solid lines and 
filled stars, and previous fits based on ~ 1 /fb of LHC data are indicated by dashed lines and open stars. 
The blue lines denote 68% CL contours, and the red lines denote 95% CL contours. 



characteristic of the NUHM1, that is largely re- 
sponsible for the expansion in the region of the 
(mo,m 1( /2) plane that is favoured at the 68% CL 
in the NUHM1, compared with the CMSSM Q 

Table [2] summarizes some interesting contribu- 
tions to the global x 2 function, namely those con- 
tributing A x 2 > 1, ATLAS 5/fb jets + # T and 
the new global combination of BR(£> S — > 
measurements. We present their contributions at 
the best-fit point in the CMSSM (which has rela- 
tively high mass parameters) , at the local best-fit 
point in the low-mass 'island' in Fig. |6j at the 
best-fit point in the NUHM1, and at a high-mass 
point in the NUHM1. We see that BR(6 -> sj) 
favours lower masses in the NUHM1, whereas 
BR(_B„ — > tv t ) applies an essentially constant 
Ax 2 penalty [^J None of the preferred models 



7 The relic density may also be brought into the WMAP 
range because the \\ acquires a relatively large higgsino 
component, another possibility made possible by the vari- 
ation in fi that is possible in the NUHM1. 
8 We note that the Belle Collaboration has recently re- 
ported a new measurement of BR(B U — > tu t ) that is in 
better agreement with the Standard Model and the classes 
of supersymmetric models discussed here [561 . 



gives a good fit to (g — 2) M , since they repro- 
duce approximately the SM value, though lower 
masses are somewhat preferred, particularly in 
the NUHM1. The lower-mass fits worsen the fit to 
M w and that in the CMSSM (NUHM1) worsens 
(improves) the fit to A^(SLD), but there are only 
small changes in crP d . We note that the ATLAS 
5/fb jets + IpT constraint favours higher masses, 
whereas the new BR(£? S — > constraint 
favours the low-mass region in the CMSSM and 
(marginally) in the NUHM1, and all the models 
fit XENON100 equally well. 

The (tan/3, mi/2) planes 

As can be seen in the left panel of Fig. [7J 
where the (tan/3, mi/2) plane in the CMSSM is 
shown, the two 'islands' visible in the (wo? TO i/2) 
plane of the CMSSM (left panel of Fig. |j) corre- 
spond to different ranges of tan j3. Large values 
of tan/3 ~ 45 are favoured only at large mi/2 > 
1500 GeV, and values of tan j3 as low as ~ 10 
come within the Ax 2 = 2.30 (68% CL) region for 
m i/2 ~ 800 GeV, consistent with their associa- 
tion with funnel and stau coannihilation regions 
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Model 


Data set 


Minimum 


Prob- 


m 


%/2 


A 


tan/3 






rVd.o.f. 


ability 


(GcV) 


(GcV) 


(GeV) 




CMSSM 


LHCi/fb 


31.0/23 


12% 


1120 


1870 


1220 


46 




ATLAS 5/fb (low) 


32.8/23 


8.5% 


300 


910 


1320 


16 




ATLAS 5 /fb (high) 


33.0/23 


8.0% 


1070 


1890 


1020 


45 


NUHM1 


LHCi/fb 


28.9/22 


15% 


270 


920 


1730 


27 




ATLAS 5/fb (low) 


31.3/22 


9.1% 


240 


970 


1860 


16 




ATLAS 5 /fb (high) 


31.8/22 


8.1% 


1010 


2810 


2080 


39 



Table 1 

The best-fit points found in global CMSSM and NUHM1 fits using the ATLAS 5/fb jets + J$t con- 
straint JW, the combination of the ATLAS J^, CDF CMS J2^ and LHCb BR(B S -> 
/i + /Lt _ ) \251] constraints and the updated values of My/ and m t , compared with those found previ- 
ously in global fits based on the LHC^ift, data set. Ln both cases, we include a measurement of 
M h = 125 ± 1.0 ± 1.5 GeV and the new XENON100 constraint J2T$. In the case of the CMSSM, we 
list the parameters of the best-fit points in both the low- and high-mass 'islands ' in Fig. [6| and we quote 
results for a high-mass NUHM1 point as well as the low-mass best-fit point in this model. We note that the 
overall likelihood function is quite flat in bot the CMSSM and the NUHM1, so that the precise locations 
of the best-fit points are not very significant, and we do not quote uncertainties. For completeness, we 
note that in the best NUHM1 fit m 2 H = m 2 H = m 2 Hd = -6.5 x 10 6 GeV 2 , compared with -5.5 x 10 6 GeV 2 
previously. 



respectively. We see again in the right panel of 
Fig. [7|that the region of the (tan/3, mi/2) plane 
that is allowed in the NUHM1 at the 68% CL 
is connected, and there is no strong separation 
between the low- and high-mass regions. As in 
Fig. [6| the region of the (tan/3, m 1( / 2 ) plane al- 
lowedat the 95% CL is also enlarged in the 
NUHM1, compared to the CMSSM. The region at 
large and relatively small tan /3 corresponds 
to the grouping of points at low mo and large 
m 1 /2 already noted in the right panel of Fig. [6] 

The (Ma, tan /3) planes 

We see in the left panel of Fig. [8j where the 
(M A , tan/9) plane in the CMSSM is shown, that 
relatively large values of Ma are favoured in this 
model, Ma > 1200 GeV at the 68% CL and 
> 1000 GeV at the 95% CL. As is seen by com- 
paring the solid contours corresponding to A% 2 
contours in this fit with the dotted lines corre- 
sponding to the analysis in [2], we see substan- 



tial retreats at small Ma and large tan/3 that 
are due to the updated BR(_B S — > n + fi~) con- 
straint. In the NUHM1, smaller values of Ma are 
allowed at low tan /3, as low as ~ 500 GeV at the 
68% CL and - 350 GeV at the 95% CL. On the 
other hand, a substantial region of the NUHM1 
(Ma, tan /3) plane at small Ma and large tan /3 is 
also excluded by the new BR(£> S — » /i + /i~) con- 
straint [3 

7. Relative Impacts of ATLAS 5 /fb, 
BR(B S // + £t~) and XENON100 

In this Section we analyze in more detail the 
different impacts of the new ATLAS 5/fb jets 
+ LpT-, BR(B S -> pL+n~) and XENON100 con- 
straints on the CMSSM and NUHM1. 

9 We see in both panels of Fig. [5] 'archipelagos' at large 
Ma that might evolve with more sampling into connected 
regions allowed at the 95% CL. In the NUHM1, this region 
is mainly populated by the points with small mo and large 
m l/2 visible in the corresponding panel of Fig.lfjl 
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Global 


33.0 


32.8 


31.8 


31.3 


■D-nEXP/SM 


1.15 


1.19 


0.94 


0.18 


■D-nEXP/SM 


1.10 


1.03 


1.04 


1.08 




9.69 


8.48 


10.47 


7.82 


M w [GeV] 


0.10 


1.50 


0.24 


1.54 


Re 


0.95 


1.09 


1.09 


1.12 


Aft (6) 


8.16 


6.64 


5.68 


6.43 


A<(SLD) 


2.49 


3.51 


4.36 


3.68 


o 

"had 


2.58 


2.50 


2.55 


2.50 


ATLAS 5/fb jets + # T 


0.09 


1.73 


0.02 


1.18 


BR(B S -> 


2.52 


1.22 


1.59 


1.70 


XENON100 


0.13 


0.12 


0.14 


0.13 



Table 2 

Summary of the contributions of the most important observables to the global x 2 function at the best-fit 
high- and low-mass points in the CMSSM and NUHM1 (those with Ax 2 > I) , and of the main updated 
observables ATLAS 5/fb jets +$ T , BR(B S -> /U + aO and XENON100. 
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Figure 7. T7ie (tan/3, mi/2) planes in the CMSSM (left panel) and the NUHM1 (right panel) including 
the ATLAS 5/fb jets + $ T constraint JBj, a combination of the ATLAS Jiffi, CDF fl^, CMS f!T 
and LHCb 24 1 constraints on BR(_B S 



H+fj,- 
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In both cases, we include a measurement of 
M h = 125 ± 1.0 ± 1.5 GeV and the new XENON100 constraint J27$. The results of the current fits are 
indicated by solid lines and filled stars, and previous fits based on ~ l/fb of LUC data are indicated by 



dashed lines and open stars, 
contours. 



The blue lines denote 68% CL contours, and the red lines denote 95% CL 
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Figure 8. The (Ma, tan /3) planes in the CMSSM (left panel) and the NUHM1 (right panel) including the 
ATLAS 5/fb jets + I$t constraint, a combination of the ATLAS, CDF, CMS and LHCb constraints on 
BR(£> S — » In both cases, we include a measurement of M^ — 125 ± 1.0 ± 1.5 GcV and the new 

XENON100 constraint '211. The results of the current fits are indicated by solid lines and filled stars, 
and previous fits based on ~ 1 /fb of LHC data are indicated by dashed lines and open stars. The blue 
lines denote 68% CL contours, and the red lines denote 95% CL contours. 



Fig. [9] displays the impacts in the {mo,m 1 / 2 ) 
planes of the CMSSM (left panels) and the 
NUHM1 (right panels) of dropping the ATLAS 
5/fb jets + Lf,T constraint (top panels), the new 
BR(£> S — > /i + /i~) constraint (middle panels), or 
the new XENON100 constraint (bottom panels). 
In each panel, we compare the 68 and 95% CL 
contours in the full global fit (dashed lines) with 
a fit dropping one of the constraints (solid lines) . 
The corresponding best fit points are shown as 
open and closed stars. 

In the CMSSM case, very little effect is seen 
when the ATLAS 5/fb jets + Ipr constraint is 
dropped. On the other hand, when BK(B S — > 
/i + is dropped there are significant extensions 
of the 68 and 95% CL regions: the two 68% CL 
'islands' at low mo merge and extend to larger m , 
and the 95% band extending to large mo broad- 
ens to a much larger range of m 1 / 2 - When the 
new XENON100 constraint is dropped, the main 
effect is the appearance of the focus-point region 
at large mo and small m^. 

In the NUHM1 case, dropping the ATLAS 5/fb 
jets + $t constraint again has little effect on the 



contours, except at small mo and m 1( / 2 - Dropping 
the BR(£? S — > /i + /i~) constraint has more effect 
on the 68 and 95% CL contours in the NUHM1. 
The effect of dropping the new XENON100 con- 
straint is qualitatively similar to that of dropping 
the new BR(_B S — > fi + fi~) constraint, except that 
it allows a larger region at small mi/2 and large 
m . 



Similarly, Fig. 10 displays the impacts in the 
(Ma, tan/?) planes of the CMSSM (left) panels 
and the NUHM1 (right panels) of dropping the 
ATLAS 5/fb jets + $t constraint (top panels), 
the new BR(B S — > /i + /i~) constraint (middle 
panels) or the new XENON100 constraint (bot- 
tom panels). As in the ^0,171^2) planes, we see 
relatively little impact in cither the CMSSM or 
the NUHM1 from dropping the ATLAS 5/fb jets 
+ IpT constraint, and a much bigger effect from 
dropping the new RR(B S — > ^ + /i~) constraint, 
particularly at small Ma and large tan (3. In both 
cases, the BR(B S — > /i + /i~) constraint reduces 
drastically the value of tan (5 at the best-fit point 
from ~ 46(27) in the CMSSM (NUHM1) to - 16. 
On the other hand, most of the impact of drop- 
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Figure 9. The (mo,mi/2) planes in the CMSSM (left panels) and the NUHM1 (right panels) displaying 
the effects of dropping the ATLAS 5/fb jets + Ifix constraint (top panels), the new BR(£? S — > fi + fi~) 
constraint (middle panels) and the new XENON100 constraint (bottom panels). In each case, the full 
global fit is represented by an open green star and dashed blue and red lines for the 68 and 95% CL 
contours, whilst the fits to the incomplete data sets are represented by closed stars and solid contours. 
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ping the new XENON100 constraint is to expand 
the preferred region at large Ma and tb < 40. 

We analyse the various x 2 contributions and 
best-fit points in more detail in Table [3j Here we 
display the principal changes in the best-fit pa- 
rameters from the new global fit (second column) 
when dropping the ATLASs/fb constraint (third 
column), dropping the new BR(_B S — > /i + /i~) 
constraint (fourth column), dropping the new 
XENON100 constraint (fifth column), and drop- 
ping all three new constraints, i.e., the analysis 
of [2] (last column). We see that it is BR(£> S — > 
H + fj,~) that, among the new constraints, has the 
most significant impact on the best-fit parame- 
ters of the NUHM1, principally tan/3. In the 
lower rows we show the \ 2 contributions of the 
direct SUSY searches, denoted as A% 2 for ^t, 
BR(£? S — > fi + and (g — 2)^. One can clearly 
observe that the contribution of (g—2)^ to the to- 
tal x 2 function increases when BR(B S — > fi + 
and the other new constraints are incorporated, 
demonstrating once again that the new data are 
pushing simple supersymmetry-breaking models 
such as the NUHM1 (and CMSSM) out of the 
'sweet spot' where the experimental measurement 
of (g — 2) M could easily be explained. 

8. Predictions for Physical Observables 

In this section we present one-dimensional x 2 
likelihood functions for sleeted physical observ- 
ables, based on the above surveys of the CMSSM 
and NUHM1 parameter spaces, comparing results 
based on the LHCi/fb and LHCs/fb data sets. 



disfavoured by Ax 2 i5 2, though the envelope of 
the x 2 function could only be reduced by further 
sampling. In the case of the NUHM1, a much 
more weakly bimodal distribution is visible. In 
the CMSSM, values of m s between ~ 1500 GeV 
and ~ 4800 GeV are allowed at the Ax 2 < 2 
level, and in the NUHM1 values over 6000 GeV 
are allowed at this level. Once again, the ability 
to move to large masses (in this case the gluino 
mass or TO1/2) can be traced to the freedom in 
the NUHM1 of adjusting Ma ~ 2m x so that the 
relic density is at the WMAP value. At such large 
masses, there is essentially no further price paid 
by g^ — 2 or other low energy observables as we 
approach the decoupling limit. 

The Squark Mass 

Fig. 12 displays the x 2 functions for the average 
mass of the supersymmetric partners of the five 
lightest right-handed quarks, mj, in the CMSSM 
(left panel) and the NUHM1 (right panel), for 



fits to the LHC 



5/fbi 



BR(£? S 



fi+fi 



and new 



XENON100 data (solid lines) _and the LHC 1/fb 

In the case of 
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data (dashed lines), as in Fig. 
the CMSSM, we see local minima of x 2 at ~ 
1800 and 3500 GeV, separated again by a range 
with Ax 2 ^ 2. At this level of A% 2 , any value of 
rriq between ~ 1200 GeV and ~ 4200 GeV is al- 
lowed. In the case of the NUHM1, there is a sim- 
ilar structure, but it is less pronounced and the 
intermediate range has Ax 2 ~ 1. Only squark 
mass values above ~ 5.4 TeV receive a substan- 
tial x 2 penalty, which is due to large ttiq being 
disfavoured in this model. 



The Gluino Mass 
Fig 
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displays the \ functions for 



the CMSSM (left panel) and the NUHM1 (right 
panel), for fits to the LHC 5 /fb data including 
the new BR(B S -> and XENON100 lim- 

its (solid lines) and to the LHCi/fb data includ- 
ing the BR(£> S — » /i + /i~) limits then available 
(dashed lines). In the case of the CMSSM, we 
see a bimodal distribution, with local minima at 
Trig ~ 2000 GeV and ~ 4000 GeV, corresponding 
to the best fits in the low- and high-mass islands 
discussed above, whose values of x 2 are almost 
indistinguishable. Intermediate values of m,g are 



The Stop Mass 

Fig. 13 displays the x 2 functions for the lighter 
stop squark mass mj in the CMSSM (left panel) 
and the NUHM1 (right panel), for fits to the 
LHC 5/fb , BR(B S -)• and new XENON100 

data set (solid lines) and the LHCi/ft, data 

We see that the 
< 1100 GeV, 
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(dashed lines), as in Fig 
CMSSM x 2 is minimized for m 
with a secondary minimum for mj ~ 2900 GeV, 
whereas in the NUHM1 x 2 is minimized for ~ 
1700 GeV with values above around 4.5 TeV be- 
ing disfavoured. In both models, values of < 
rriq are expected. 
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Figure 10. The (Ma, tan (3) planes in the C'MSSM (left panels) and the NUHM1 (right panels) displaying 
the effects of dropping the ATLAS 5/fb jets + Ifix constraint (top panels), the new BR(£? S — > fi + fi~) 
constraint (middle panels) or the new XENON100 constraint (bottom panels). In each case, the full 
global fit is represented by an open green star and dashed blue and red lines for the 68 and 95% CL 
contours, whilst the fits to the incomplete data sets are represented by closed stars and solid contours. 
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Quantity 


.bull lit 


w/o AlLAb 5 /fb 


w/o B s — > fi + /i 


w/o XrjJNUlNlUO 


T 1 1 ' £+- 

LitlKjl/fb nt 


Total x 2 


31.3 


30.1 


28.7 


31.3 


28.9 


m (GeV) 


240 


240 


270 


240 


270 


m 1/2 (GeV) 


970 


870 


920 


970 


920 


A (GeV) 


1860 


1860 


1730 


1860 


1730 


tan/3 


16 


15 


27 


16 


27 


m% (GeV 2 ) 


-6.5 x 10 6 


-5.6 x 10 6 


-5.5 x 10 6 


-6.5 x 10 6 


-5.5 x 10 6 


A X 2 # T 


1.18 


1.06* 


1.36 


1.18 


0.85* 


Ax 2 B s -> 


1.70 


1.12 


0.21* 


1.70 


0.21* 


Ax 2 XENON100 


0.13 


0.13 


0.13 


0.14* 


0.14* 


AX 2 (3 - 2) M 


7.82 


7.41 


5.99 


7.82 


5.99 



Table 3 



Principal changes in the best-fit parameters in the NUHM1 from the new global fit (second column) when 
dropping the ATLAS^/f^ constraint (third column), dropping the new BR(_B S — > /i + /i~) constraint (fourth 
column), dropping the new XENON100 constraint (fifth column) and dropping all three new constraints, 
i.e., the analysis of J3j (last column). The bottom four lines display the contributions to \ 2 °f ^ e 
JfiT: BR(£? S — > /i + /i~), XENON100 and {g — 2) M constraints. The numbers marked by * are for the 
implementations of these constraints in Wl. 




Figure 11. The x 2 likelihoods in the CMSSM (left panel) and NUHM1 (right panel) as functions of nig 
based on global fits to the LHC^/f^, BR(_B S — > /i + /i~) and new XENON100 data set (solid lines), and to 
the LHCx/fo data set (dashed lines). 
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Figure 12. The x 2 likelihoods in the CMSSM (left panel) and NUHM1 (right panel) as functions of mq 
based on global fits to the LHC^/^, BR(B S — > fi + /i~) and new XENON100 data set (solid lines), and to 
the LHCi/ib data set (dashed lines). 




Figure 13. The \ 2 likelihoods in the CMSSM (left panel) and NUHM1 (right panel) as functions of 
based on global fits to the LHCs/ib, BR(£> S — > and new XENON100 data set (solid lines), and to 

the LHCi/fo data set (dashed lines). 



The Stau Mass 
Likewise, Fig. 



14 



displays the \ 2 functions for 
m fl in the CMSSM (left panel) and the NUHM1 
(right panel), for fits to the LHC 5 /fb, BR(B S — > 
and new XENON100 data set (solid lines) 



and the LHCi/fb data (dashed lines), as in Fig.fTTj 
In this case, in the CMSSM % 2 IS minimized 
800 GeV, whereas in the 



for mfj <~ 400 and 



NUHM1 x is minimized for m fl - 700 GeV with 
an upper limit around 2.5 TeV. 
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4000 5000 



m Tl [GeV| m Ti [Gey| 



Figure 14. The x 2 likelihoods in the CMSSM (left panel) and NUHM1 (right panel) as functions of 
based on global fits to the LHC^/f^, BR(_B S — > and new XENON100 data set (solid lines), and to 

the L-ffCi/fb data set (dashed lines). 



The Higgs Mass 

In Fig. 15 we show as solid lines the % 2 func- 
tions for M h in the CMSSM (left panel) and the 
NUHM1 (right panel), for fits to the LHC 5/fb , 
BR(B S -> n+pr) and new XENON100 data set, 
but omitting all the direct constraints on the 
Higgs mass from the LHC, LEP and Tevatron. 
Thus, Fig. [15] displays the predictions for Mh in 
the CMSSM and NUHM1 on the basis of the 
other constraints. As in previous publications 
(see [T][2l[T4] and references therein), the plots 



in Fig. [15] include a red band with a width of 
1.5 GeV to represent the theoretical uncertainty. 
These plots show directly the extent to which 
the indirect prediction of Mh agrees with the 
(provisional) direct measurement. Also shown in 
Fig.[l5]is the LEP exclusion for a SM Higgs boson 
(shaded yellow), and the range indicated by the 
LHC discovery (shaded green). We see that the 
LHC measurement of is somewhat larger than 
would be indicated by our best fit based on other 
data, but not disastrously so, with A\ 2 ~ 1.5 in 
the CMSSM and - 0.5 in the NUHM1. 

BR(B S -> fi + n~) 

Fig. 16 displays the x 2 functions for BK(B S — > 
M+A* - ) in the CMSSM (left panel) and the 



NUHM1 (right panel), for fits to the LHC 5/f b, 
BR(B S -+ n+yr) and new XENON100 data set 
(solid lines) and the LHCi/ft, data (dashed lines), 
as in Fig. [TT] Comparing the solid and dot- 
ted lines, we see the strong impact of the recent 
BR{B S -> m + aO update. In both the CMSSM 
and the NUHM1, values of BR(B S -> fi + fj,~) 
above the SM value are preferred over values be- 
low the SM value. We see that, in either case, a 
measurement of BK(B S —> /i + /i~) at the level of 
30%, which may soon be possible, would greatly 
restrict the preferred ranges of both the CMSSM 
and NUHM1 parameter spaces. 

ct si 

P 

Fig. [17] displays the best-fit points and the 
68% and 95% CL contours (red and blue, respec- 
tively) in the (m^o , crjj 1 ) planes in the CMSSM 
(left panel) and and the NUHM1 (right panel). 
The solid lines are for the global fit to the LHCs/ft, 
BR(B S -> n+fjT) and new XENON100 data set, 
and the dashed lines for the previous LHCiyfb fit, 
and the solid (open) green stars mark the corre- 
sponding best-fit pointsF^l We see that the region 



u We display results for m~o < 1 TeV only, because re- 
sults from XENON100 are not published for larger masses. 
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Figure 15. The one- dimensional Ax 2 functions for Mh in the CMSSM (left) and the NUHM1 (right). 
The solid lines are for fits including all the available data including (g — 2) M but excluding the direct 
LEP [61, 62 1 , Tevatron [?] and LHC \36\[3'?$ constraints on Mh, with a red band indicating the estimated 
theoretical uncertainty in the calculation of Mh of ~ 1.5 GeV. The shaded ranges are those excluded by 
LEP searches for a SM Higgs boson (yellow) and indicated by the LHC discovery (green). 
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Figure 16. The \ 2 likelihoods in the CMSSM (left panel) and NUHM1 (right panel) as functions of 
BR(S S — > fi + [i~) based on global fits to the LHC^/^, BR(_B S — > and new XENON100 data set 

(solid lines), and to the LHCi/^ data set (dashed lines). 



of the (m^o,^ 1 ) favoured in the CMSSM is now were previously allowed at the 95% CL. 
more restricted, in particular at small rn^o. The 
main effect of the new XENON100 constraint has 
been to remove focus points with large a? 1 that 
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m^o[GeV] m x o[GeV] 



Figure 17. The (m%o, a^ 1 ) planes in the CMSSM (left panel) and the NUHM1 (right panel). The Ax 2 = 
2.30(5.99) contours, corresponding to the 68(95)% CL are coloured red (blue). The solid (dashed) lines 
are for global fits to the LHC^/^, BR(_B S — > an d new XENON 100 (LHC±/fo) data, and the 

corresponding best-fit points are indicated by solid (open) green stars. 



9. Comparison with other Analyses 

We are not aware of any global analyses in- 
corporating the new XENON100 constraint [27] . 
The two global analyses including LHC 5/fb data 
that are most similar to ours are |] and [10], 
and we take this opportunity to comment on the 
relations and differences between those and the 
present work. 

As in this work, Ref. [9] incorporates the con- 
straints imposed by the LHC 5 /fb data using an 
implementation of Delphes that is reported to re- 
produce the 95% CL limits published previously 
by ATLAS using 165/pb and 1/fb of data and 
shown to reproduce approximately the 95% CL 
limits for tan/3 — 10 and Aq found in the 
ATLAS analysis used here and the CMS razor 
analyss with 5/fb of data 13 . However, the 
ATLAS 5/fb jets + Jf*T analysis has apparently 
not been modelled in detail, and full informa- 
tion is not available on its validation for other 
values of the CMSSM parameters, nor on its ex- 
tension to the NUHM1. In contrast to previous 
papers by the same collaboration that used the 
Master Code, Ref. [9] uses only a reduced set of 
precision electrowcak and flavour observables, a 
topic we comment on below. In contrast, Ref. [9] 



uses AstroFit, which provides input from both 
direct and indirect dark matter searches, whereas 
we use only the direct XENON100 search. How- 
ever, Ref. [9j finds that the indirect dark mat- 
ter searches have negligible impact at the present 
time. When the constraint M h = 126±2±3 GcV 
is imposed in [9], the results for the CMSSM are 
relatively similar to ours, whereas the results for 
the NUHM1 are only qualitatively similar. 

Ref. [To] implements the CMS razor analysis 
using 5/fb of data, and provides more details of its 
efficiency and likelihood maps, and also compares 
with the 95% CL limit obtained in the ATLAS 
5/fb jets + $t analysis. This paper favours val- 
ues of the CMSSM mass parameters (7710,^1/2) 
that are considerably larger than in our analysis 
and in [9] . However, Ref. [10] does favour the stau 
coannihilation strip and the rapid-annihilation 
funnel regions of the CMSSM, disfavouring the 
focus-point region that had been advocated in 
previous fits by the same group. Ref. [lO] does not 
discuss the NUHM1, but does consider the possi- 
bilities that \i < and of dropping the (g — 2) M 
constraint. 
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10. Conclusions 

We have presented in this paper new global fits 
to the CMSSM and NUHM1 based on the 7 TeV 
ATLAS 5/fb jets constraint and the latest 

global combination of limits on BR(_B S — > fj, + n~), 
the new XENON100 constraint, a Higgs mass 
measurement Mh ~ 125 GeV, and minor updates 
for other observables. Before the advent of the 
LHC, low-energy data, in particular (g — 2) M , gave 
hope that the scales of supersymmetry breaking 
in these models might be within reach of the LHC 
within its first year of operation at 7 TeV in the 
centre of mass. This has not been the case, and 
the p- values of the CMSSM and NUHM1 have de- 
clined to 8.5 and 9.1%, respectively, prompting a 
re-evaluation of their status. 

Two possibilities remain open at this point: ei- 
ther supersymmetry is relatively light but man- 
ifests itself differently from the CMSSM or 
NUHM1, or the supersymmetry-breaking mass 
scale is heavier, possibly even beyond the reach 
of the LHC. In the first case, several possibili- 
ties suggest themselves (see, e.g. 



Ref. 63 ): other 



models of supersymmetry breaking, light third- 
generation squarks, J"-SU(5) (D], the pMSSM, 
R- violating models, etc., and a discussion of these 
possibilities lies beyond the scope of this work. As 
for the second case, the global fits presented here 
in the CMSSM and NUHM1 give limited guid- 
ance on the possible supersymmetric mass scale. 

Perhaps surprisingly, the ATLAS 5/fb jets + 
JpT data are not much more problematic for the 
CMSSM and the NUHM1 than were the LHC 1/fb 
data in combination with the Higgs mass con- 
straint. The latter two had already constrained 
the CMSSM and NUHM1 parameters to lie be- 
yond the range where they could fully 'solve' the 
(g — 2) M problem, and the contribution of the AT- 
LAS 5/fb jets + IpT constraint to the global \ 2 
function is not large, even at the local best-fit 
point within the low-mass CMSSM 'island'. A 
new feature of this analysis is the impact of the 
latest BR(£? S — > constraint, and we stress 

the importance of the constraint from the search 
for spin-independent dark matter scattering by 
the XENON100 Collaboration [27]. 

The most important upper limit on the 



supersymmetry-breaking mass scales in the 
CMSSM and the NUHM1 is now provided by the 
cold dark matter density constraint, and indica- 
tions on the possible sparticle masses within the 
range allowed by this constraint are now quite 
weak. The (g — 2) M constraint still offers some 
preference for low sparticle masses, but this is 
largely counterbalanced by other constraints such 
as BR(£? S — > As we have seen, any 

value of m~g € (1500, 5000) GeV is allowed at the 
Ax 2 < 2 level. The lower part of this range may 
be explored with LHC operating at 8 TeV, and 
the middle part with the LHC at 13 or 14 TeV, 
but the upper part of the allowed range of nig 
would lie beyond the reach of the LHC. 

On the other hand, one must keep in mind that 
the increasing tension within the CMSSM and 
NUHM1 between the low-energy data, mainly 
(g - 2) p and the LHC5/fb and BR(B S -> ^ + A*~) 
data on the other side, has resulted over time in 
a monotonously-increasing x 2 and a correspond- 
ingly lower fit probability. Consequently, looking 
beyond the CMSSM and NUHM1, to models with 
a different connection between the coloured and 
uncoloured sector, not only seems timely now, but 
mandatory. We expect that some model with 
such a different connection will give a substan- 
tially better fit to all existing experimental data 
and present more favorable expectations for sub- 
sequent LHC runs or a future LC. 
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